Operator fidelity susceptibility: an indicator of quantum criticality 
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We introduce an operator fidelity and propose to use its susceptibility for characterizing the 
sensitivity of quantum systems to perturbations. Two typical models are addressed: one is the 
transverse Ising model exhibiting a quantum phase transition, and the other is the one dimensional 
Heisenberg spin chain with next-nearest-neighbor interactions, which has the degeneracy. It is 
revealed that the operator fidelity susceptibility is a good indicator of quantum criticality regardless 
of the system degeneracy. 
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Introduction — There are two important concepts, en- 
tanglement and fidelity in quantum information the- 
ory These two concepts are closely related to each 
other. For instance, fidelity, which was first proposed 
as a tool for describing the stability of a quantum sys- 
tem to perturbations [3L may be used to characterize 
quantum entanglement Q . Notably, fidelity has recently 
been used to characterize quantum phase transitions 
(QPTs) 0, H, H, 0| ■ On the other hand, entanglement 
has also been employed to be an indicator of QPTs in 
many correlated quantum systems [1, d, [l(| E3 • 

How to characterize the stability of a quantum sys- 
tem to perturbations is an important issue as there is 
no quantum counterpart of the classical Lyapunov ex- 
ponent. The Loshmidt echo [l2| has been adopted as 
a measure of the system stability against perturbations, 
which is introduced as follows. Let operators Uq and Ui 
denote the time evolutions of Hamiltonians Hq and H\ , 
where Hi is slightly different from Hamiltonian Hq with 
Hi — Hq = eV as a small perturbation. In this case, the 
operator U c = UqIIi is referred to as the echo operator, 
and the absolute value of its expectation over a specific 
state is defined as the Loshmidt echo 



L m = mulUi\vb)\. 



(1) 



This is just the fidelity amplitude. Obviously, it is state- 
dependent, i.e., one has to choose an initial state (artifi- 
cially in many cases) to evaluate its response to pertur- 
bations. This scenario to characterize QPT has a serious 
limitation, e.g., it can hardly be applied to a degenerated 
ground state, which has been a great challenge for a long 
time. 

In this Letter, mainly motivated by the above chal- 
lenge, we introduce a new kind of fidelity measure, called 
operator fidelity, and propose for the first time to use 
its susceptibility for characterizing the stability of quan- 
tum systems to perturbations. A distinct and significant 
merit lies in that it is state-independent and, in partic- 



ular, is able to characterize the quantum criticality re- 
gardless of degeneracy. To illustrate the feasibility and 
reliability as well as the merit of the introduced opera- 
tor fidelity susceptibility, here we employ it to investigate 
two typical QPT systems: the quantum Ising model and 
the Heisenberg model with next-nearest-neighbor inter- 
actions. Indeed, this fidelity susceptibility is able to serve 
as an indicator of QPT. In addition, for comparison, we 
also consider the mixed state fidelity susceptibility to ad- 
dress the quantum criticality with the ground-state de- 
generacy. 

We begin with the definition of operator fidelity. Let 
Ti be a c£-dimensional Hilbert space. All linear operators 
on Ti are represented by d x d matrices and thus their 
own may be considered to be vectors in an expanded 
eP-dimensional Hilbert space Tins- The inner product 
Tins is defined as the Hilbert-Schmidt product, i.e., for 
operators A and B, (A\B) — Ti(A^B). In this sense, 
any linear operators on Ti. can be considered as a state 
on Tins- Thus, the fidelity of two states can naturally 
be generalized to the operator level. For two unitary 
evolution operators Uq and Ui on Ti, the fidelity between 
them is defined as 

F 2 = ±\Tr(tf Ui)\ 2 = mutUi)] 2 , (2) 

where the averaged tracing operation is defined as 
Tr () =Tr() jd. It is notable that one may obtain the av- 
eraged Loshmidt echo [IH-pjl after averaging Lim over 
all states on Ti with a Haar measure, and the averaged 
Loschmidt echo and the operator fidelity are essentially 
equivalent. Remarkably, the operator fidelity involves 
not only the ground state, but also all eigenstates of the 
system. It quantifies the difference between two unitary 
operators, and is a conserved quantity under local oper- 
ation in Hilbert space Tins- 

We can rewrite the echo operator U e as [13] 



U e = 1-ie 



Vi(h)dti 



2 



-e 2 / / V / (ii)y / (t 2 )dt 1 di 2 + 0(e 3 ), (3) 
Jo Jo 

where Vr(t) = exp(ii/oi)^(*) ex P( — i-Hio*) is the pertur- 
bation operator in the interaction picture. After tracing, 
we have 



Tr(t/ C ) = 1 - ieTr[W(t)\ - -Tr[W(t) 2 ] + 0(e 3 ), (4) 
where W(t) = J* Vi(t')dt'. Then from Eq.(2), we obtain 

F 2 = 1 - e 2 [T7(W/(i) 2 ) - T^ 2 (W(t))} + 0(e 4 ). (5) 

To evaluate the above operator fidelity, one has to choose 
a small parameter artificially, which is e-dependent. To 
avoid this artifact, we can also introduce a so-called fi- 
delity susceptibility [J, [l8[ , which is given by 



Xr- 



lim 

E^O 



1 - F 



i{Tr [W{tf] 



Tr [W{t))}. (6) 



Remarkably, the above simple formula possesses a dis- 
tinct computational advantage that enables one to cal- 
culate straightforwardly the fidelity susceptibility from 
W(t), which can also be evaluated readily or at least nu- 
merically for more complicated systems. On the other 
hand, generally speaking, a quantity/measure suscepti- 
bility responds to the relevant perturbations more sensi- 
tively than the quantity/measure itself does, so we be- 
lieve that it could capture a drastic change feature of 
the quantum evolution (versus the relevant parameter) 
around a critical point. We below explore the intriguing 
relationship between the operator fidelity susceptibility 
and the QPT in two typical systems, with one having 
degeneracy. 

Quantum phase transition — The first system we con- 
sider is an Ising spin chain subject to a transverse mag- 
netic field, whose Hamiltonian reads 



Hn = 



M 

£ 

l=-M 



(7) 



where A characterizes the strength of the transverse field, 
erf (a = x, y, z) are the Pauli operators defined on the l- 
th site, and the total number of spins in the Ising chain is 
N = 2M +1. The perturbation operator is given by eV — 
e^2lf._M erf /2. There are two competing terms in the 
Hamiltonian, i.e, the Ising interaction and the transverse 
field term. 

The Hamiltonian can be diagonalized by combining 
Jordan- Wigner transformation and Fourier transforma- 
tion to the momentum space, i.e., 



(Slfccr/cz) e" 



fe>0 



l-£)*a. (8) 



where we have used the following pseudospin 
operators a ka (a = x, y, z): cr kx 



d-kdk, (k — 1, 2, ...M) , o~k y = ~id\d}_ k + id_ k d k ,cr kz = 

d\d k + d)_ k d- k - 1,00s = 2d d - 1. Operators 

d\,dk{k = 1,2, ...M} denote the fermionic creation and 
annihilation operators in the momentum space. Here, 



-A+2cos 



/2vrfc 



4 sin 



2irk 



Ok = arcsin 



(9) 



Then the time evolution operator is derived as (with 

h = i) 

fc>0 

The unitary operator U\{t) for Hamiltonian H± = Hq + 
eV can be obtained by just replacing A with A + e in the 
above equation. 

At this stage, from Uo(t) and Ui(t) given above, we 
are able to obtain W(t) as, 



W(t) = ^[(T fez (tcos 



fe>0 



^9 fe + ^sin2ffi fe 



+a ky cos 9 k sin 6 k [t - — — sin 2tfl k 
V 2il k 

+<r kx smO k —— (cos2ffi fc - 1)] + ^-a 0z (U) 
2\l k 2 

Consequently, the fidelity susceptibility is derived exactly 




cos Qi 



fe>0 



-J2 s[n2 ( n kt) sin 2 O k /n 



(12) 



fe>0 



Note that the first term in Eq. (|12|) . which is proportional 
to the square of time t, plays a dominant role when t is 
large. 

In the transverse Ising model, two phases are sepa- 
rated by the quantum phase transition point A = 2. The 
singular behavior of QPT at the transition point reflects 
the sensitivity of ground state to perturbations. At this 
stage, we numerically look into the behaviors of the oper- 
ator fidelity susceptibility and its partial derivative with 
respect to A at a finite time t — 100 (the natural units are 
used here). As shown in Fig. 1 for different system sizes, 
the transition point is unambiguously signatured: it is 
clearly seen that the fidelity susceptibility and its partial 
derivative are nearly unchanged when increasing A from 
to 2; the derivative increases sharply at the transition 
point and the derivative peak is higher when the system 
size becomes larger. 

The behavior of the operator fidelity susceptibility at 
the QPT point is different from that of the ground-state 
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FIG. 1: Partial derivative of the fidelity susceptibility versus 
the parameter A for different system sizes TV = 2 10 ,2 11 ,2 12 . 
The insert plots the fidelity susceptibility xf versus the pa- 
rameter A. The time t = 100. 

fidelity The ground-state fidelity susceptibility dis- 
plays a sudden increase at the QPT point, reflecting 
the drastic change of ground state of the system when 
the QPT occurs, while the operator fidelity susceptibil- 
ity drops to a certain value continuously at the QPT 
point (and is unchanged below the point) since it involves 
the all eigenstates and characterizes the sensitivity of 
the whole system to perturbations in the time evolution. 
However, on the other hand, its partial derivative changes 
discontinuously at the QPT point and is much more sen- 
sitive to perturbations, being able to single out the QPT 
point unambiguously. 

Heiserberg model with next-nearest-neighbor 
interactions — For the fidelity scenario developed 
previously for QPTs, only pure ground states can be 
addressed, without taking into account the degeneracy; 
while it is the case for some quantum systems. As 
seen above, the operator fidelity approach has an 
advantage that the degeneracy is not necessary to 
be considered explicitly. To contrast our approach 
with the state fidelity approach, we below address a 
model with the ground-state degeneracy. The Hamil- 
tonian of one-dimensional Heisenberg system with 
next-nearest-neighbor interaction reads 

N 

H = (JiSj • s i+ i + J 2 Sj • s l+2 ) , (13) 

i 

where the Sj denotes the spin- 1/2 operator at the ith 
site, TV is the total number of sites, J\ and J 2 are the 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
exchange couplings. As usual, we choose the periodic 
boundary condition and set J% = 1 for convenience. The 
perturbation operator is given by eV = e^jSj • s^ +2 . 
Note that no exact analytical results are available for 
this model (fTB")) except the special case of J 2 = and 
J 2 = 1/2. 

It is well known that the point J 2 = 1/2 corresponds 
to the Majumdar- Ghosh model where the ground state 
is the products of dimers, leading to a gaped phase [2l| . 




FIG. 2: (a) The operator fidelity susceptibility versus J2 for 
different system sizes TV = 7,9,11. (b) The mixed state fi- 
delity susceptibility versus J2. 

Chen et al [13] studied the ground-state fidelity and first- 
excited-state fidelity of this system with even number of 
sites. Here we focus on the odd number of sites as the 
fourfold degenerate energy level structure is present in 
this case. 

We first diagonalize the Hamiltonians numerically, and 
then calculate the operator fidelity susceptibility versus 
the NNN coupling J 2 for TV = 7, 9, 11, as plotted in Fig. 
2(a). At a finite time t = 100, the susceptibility \ F de- 
cays to a minimum value near the critical point J 2 = 0.5. 
With the size increasing, the minimum point is closer to 
the critical point. It is expected that the curve around 
the minimum point would become sharper and sharper 
when the size increases, leading to a discontinuity in its 
partial derivative with respect to J 2 at the critical point 
in the thermodynamic limit, as in the case of the trans- 
verse Ising system. We indeed note from the energy spec- 
trum that the ground energy level and the excited energy 
level crosses near the point J 2 = 0.5. In this sense, the 
operator-fidelity susceptibility (or its partial derivative) 
is also able to capture the level crossing feature in the sys- 
tem and thus to indicate the critical point, overcoming 
the subtle problem induced by the degeneracy. 

On the other hand, at least for comparison, it is also 
interesting to consider an alternative approach to address 
degenerate cases by making use of the mixed state fidelity 
given by [23| 

F (p , Pl ) = Tr (^Jp\ /2 pvp\ ,2 \ = Tr (^pj) . (14) 

Without loss of generality, it is not unreasonable to as- 
sume the mixed ground state as an equal mixture of the 
degenerate ground states, 

1 R 

Pi = ft22\^jr)(lpjr\, (15) 
r— 1 

with r = 1,2, ...R denote the degeneracy and the state 
\4'jr) denotes the jth degenerate eigenstate of the system. 
In the fidelity F (po> Pi) of this Heisenberg spin chain 
with the NNN interactions, po comes from the mixture 
of the ground states of Ho, and p\ corresponds to H\ = 
Ho + eV. 



4 



When the degeneracy of the system is explicitly ob- 
tained, we evaluate the mixed state fidelity susceptibility 
versus the coupling strength J2 by combining Eq. (6) with 
Eq. (14), as shown in Fig. 2(b). Clearly, the susceptibil- 
ity x F passes the critical point J 2 — 0.5 discontinuously. 
For larger sizes such as N — 9, 11, there exist two peaks 
since the ground energy level crossing occurs twice. With 
increasing the system size, the position of the first peak 
approaches to the critical point, and the second one is 
closer to the first one. Although, it seems that the sug- 
gested mixed state fidelity approach may also indicate 
the critical point, it should be pointed out that it is fea- 
sible only when the degeneracy of the ground states is 
explicitly known and the equal mixture of the degener- 
ated states is assumed. In addition, due to the energy 
level crossing, the degeneracy may change at the critical 
point and thus the mixed state fidelity approach may not 
be workable for all the values of the considered parame- 
ter. 

Relationship to entangling power — Finally, we would 
like to disclose an intrinsic connection between the 
present operator fidelity and the entangling power [lil] 
that was adopted to characterize the entangling capa- 
bility of a quantum evolution. The entangling power is 
essentially the mean state linear entropy at time t after 
averaging over all initial states. We consider a general 
Hamiltonian in the form 

H = I®H + \1)(1\®V = |0)<0|<g)E/o + |l><l|<S>#i, (16) 

where / is the identity operator and Hi = Hq + 
V. The time evolution operator is readily obtained as 
U(t) = |0)(0| <g> U (t) + 1 1 > < 1 1 <g> U x {t), where U k (t) = 
cxp(—iHkt){k = 0, 1}. This is a kind of the controllcd- 



U operator, for which the entangling power e „ is pro- 
portional to the operator entanglement E[U(t)] [20], i.e., 
e p [U(t)} = [d/(d+ l)] 2 E[U{t)]. The operator entangle- 
ment E[U (t)] can be straightforwardly obtained from the 
expression of U(t). Finally, we obtain 

This establishes a direct connection between the entan- 
gling power and operator fidelity. 

Summary — We have proposed an operator fidelity ap- 
proach to characterize the stability of quantum system to 
perturbations, which possesses a remarkable advantage 
that it is state-independent and is able to characterize the 
quantum criticality regardless of degeneracy. We have 
employed the approach to reveal successfully the QPT 
points in two typical systems: the quantum Ising model 
and the Heisenberg chain with next-nearest-neighbor in- 
teractions, with the latter having the degeneracy. Our 
approach is quite promising for the exploration of quan- 
tum instability including QPT and quantum chaos. 
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